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Abstract 



In this work we refine the method of and obtain a novel kind of 
functional equation determining the partition function of the elliptic 
SOS model with domain wall boundaries. This functional relation is 
provenient from the dynamical Yang-Baxter algebra and its solution is 
given in terms of multiple contour integrals. 
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1 Introduction 

Face models or Solid-on-Solid (SOS) models of statistical mechanics were introduced by 
Baxter in the process of solving the eight- vertex model with periodic boundary conditions 
[2]. The Boltzmann weights of Baxter's eight- vertex model are parameterised by elliptic 
functions and this feature is intrinsically connected with the requirement that the model 
statistical weights satisfy the Yang-Baxter equation (3j |4j . The elliptic nature of the 
eight-vertex model Boltzmann weights is naturally transported to the corresponding 
SOS model and a new continuous parameter emerges in the course of Baxter's vertex- 
face transformation [2] . We shall refer to this new parameter as dynamical parameter [2] 
and its implications for the analytic theory of the eight- vertex model have been discussed 
in [5]. Besides the emergence of this new parameter, the resulting statistical weights no 
longer satisfy the standard Yang-Baxter equation but its dynamical version introduced 
by Felder [6]-[8] as the quantised form of a modified classical Yang-Baxter equation (9 , 10 



In the same fashion as Drinfeld-Jimbo quantum groups 11 14 provide the algebraic 
structure underlying the solutions of the Yang-Baxter equation, the so called elliptic 
quantum groups introduced in [6j[7] accomodate the solutions of the dynamical Yang- 
Baxter equation. In this work we shall restrict ourselves to the SOS model built out of 
the solution of the dynamical Yang-Baxter equation associated with the elliptic quantum 
group E Ttl [s\.2\. As far as the boundary conditions are concerned, we shall consider the 
case of domain wall boundaries firstly introduced by Korepin in the context of vertex 
models 15 and subsequently extended for SOS models in (T 16 18 . 

In contrast to the case with periodic boundary conditions, the partition function of 
vertex and SOS models with domain wall boundaries can be exactly computed without 
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relying on solutions of Bethe ansatz equations. Interestingly enough, the exact solu- 
tion of the six-vertex model with domain wall boundaries fi~9l revealed that the model 



free-energy differs from the case with periodic boundary conditions 20 . This unusual 
dependence of bulk thermodynamic properties with boundary conditions has also been 
observed for the elliptic SOS model when the anisotropy parameter assumes a particular 



value 21 . For general values of the anisotropy parameter this study still poses as an 
open problem, probably due to the lack of suitable expressions for the partition function 
allowing to compute physical properties in the thermodynamic limit. In searching for 
alternative representations for this partition function, which might render the analysis 
of the thermodynamic limit feasible, we have obtained in fl] a multiple integral formula 
for the partition function of the trigonometric SOS model with domain wall boundaries. 
This case consists of a particular limit of a more general elliptic model, the limit where 
elliptic theta-functions degenerate into trigonometric functions, and here we refine and 
generalise the method of [I] for the general elliptic case. 

This paper is planned as follows. In the Section [2] we give a brief description of SOS 
models with domain wall boundaries in terms of the generators of Felder's dynamical 
Yang-Baxter algebra. The conventions employed here are basically the ones already 
discussed in jlj. In the Section|3]we demonstrate how the dynamical Yang-Baxter algebra 
can be explored in order to obtain a functional equation determining the model partition 
function. This functional equation is solved in Section [4] and concluding remarks are 
discussed in Section [5j Technical details required throughout this paper are presented in 
the Appendices. 



2 Operatorial description of the SOS model 



Partition functions of two-dimensional lattice models can be described in terms of op- 
erators representing the allowed configurations of the lattice. This feature goes back to 



Kramers and Wannier transfer matrix technique 22,23 and it has found several impor- 



tant generalisations 24 . Remarkably, when the statistical weights of the model satisfy 
the Yang-Baxter equation [3] or its dynamical counterpart [6-8], we not only have an 
operatorial description of the model but also an algebra governing its operators for any 
size of the lattice. In what follows we shall recall the conventions discussed in [l] which 



consist of an extension of the ones given in 15 for the six- vertex model 



Dynamical Yang-Baxter equation. Following [7j|8] we encode the statistical weights 
of our elliptic SOS model on a matrix 1Z G End(V <8> V) with V = C 2 . For variables 
Aj,7, 9 G C, this matrix 1Z satisfies the dynamical Yang-Baxter equation, 



k 12 {\i -x 2 ,e- ih)n lz {\ x - a 3 , 9)n 23 (\ 2 -x 3 ,e- 7 Ao 

K 23 (\ 2 - A 3 , 0)fti 3 (Ai - A 3 , 9 - 7^ 2 )^i 2 (Ai - A 2 , 



(2.1) 



where h = diag(l, —1). The Eq. (2.1) is defined in End(Vi <S> Y 2 <S> V 3 ) and the action of 
7^i 2 (A, 6 — -yh 3 ) on the basis vector v 1 <g> v 2 <S> v 3 is understood as 



[R(X, 9 - 7%i <g> v 2 ] ® v 3 



(2.2) 
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where h is a scalar denoting a particular eigenvalue of h, i.e. h(Vi = hv^. 



Definition. Let r be a complex number such that Im(r) > and write p 
that \p\ < 1. For A G C we define the elliptic function / with nome p as 



1 oo 

/(A) = - £ (-i) n -^ 



2 p("+§) 2 e -(2n+l)A 



(2.3) 



The function /(A) corresponds to the Jacobi theta-function <9i(iA,r)/2 |25j 
The equation ( |2.1 ) has been considered in |6}[7j and its solution reads 



with non-null entries 



fa+(X,9) 



V o 



a±(A,0) 
b ± {\,9) 

c±(\,9) 











b + (X,9) c + (A, 
c_(A,0) 6_(A, 




/(A + 7) 
,/(^T) 



o \ 




a_(A,0)/ 



/(A) 
/(7) 



JgrA) 



(2.4) 



(2.5) 



The algebraic structure underlying (2.4) is the elliptic quantum group i^^ffffe] and we 
have collected the properties of / required through this work in the Appendix |Bj 



Dynamical monodromy matrix. Let d{ be the operator valued parameter 

L 

k=i+l 

and consider the following ordered product of dynamical 7£-matrices, 



(2.6) 



r a (x,e) = l[n ai (x-fi l ,e l 



(2.7) 



i=l 



living in the tensor product space V a (g) Vi <E> ■ ■ ■ <E> Vl. We shall refer to T a (X,9) as 
dynami cal m onodromy matrix or simply monodromy matrix. Since the dynamical 1Z- 
matrix (2.4) satisfy the weight-zero condition [TZ a b(X, 9),h a + h b ] = 0, one can show that 



(2.7) obeys the relation 



^(Ai - X 2 ,9 - 7 H)7;(Ai, 9)T b (X 2 , 9 - jh a ) = %(X 2 , 9)T a (X h 9 - jh b )K ab (Xi - A 2 , i 



(2.8) 
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with H = y^^fc- The relation (2.8) shall be referred to as dynamical Yang-Baxter 
fc=i 

algebra and since V = C 2 , the dynamical monodromy matrix can be recasted in the 
form 

'A(X, 9) B(X, 
C{X,9) D(X, i 



T a {X,9)-. 

whose entries are then defined on Vi 



(2.9) 



Domain wall boundaries. The partition function of the elliptic SOS model with 
domain wall boundaries can be written in terms of entries of (2.9) as described in [I]. 
More precisely, the elliptic SOS model partition function Z e is given by the expected 
value 



Z e = (0\l[B(X J ,6+n)\0) 

3=1 



where 



10} 



and 



10} 



i=l 



(2.10) 



(2.11) 



In the next section we shall demonstrate how the dynamical Yang-Baxter algebra can 
be employed to produce a functional equation determining Zq. 



3 Functional relations 



The relation (2.8) encodes commutation rules for the operators A(X,6), B(X,6), C(A, 0) 



and -D(A, 9) once the structure (2.9) is considered. Out of the sixteen relations contained 



in (2.8), we will make use of only two of them in order to derive a functional equation 



describing the partition function (2.10). More precisely, the required relations are simply: 



B{X 1 ,9)B{X 2 ,9 + 1 ) 
A(X U 9 + 1 )B(X 2 ,9) 



B(X 2 ,9)B(X 1 ,9 + 1 ) 
/(A2-A1 + 7) f{0 + l) 
/(A2-AO /(0 + 2 7 ) 
/(^ + 7 -A 2 + A 1 ) /( 7 ) 



/(A 2 - Ax) f(6 + 2 7 ) 



B(A 2 ,0 + 7)A(Ai,0 + 27) 

B(X 1 ,9 + 1 )A(X 2 ,9 + 2 1 ) . 



(3.1) 



In addition to that, the weight-zero condition satisfied by (2.4) associated with the 



definition (2.7) allow us to compute the action of A(X, 9) on the states |0) and |0) defined 



in (2.11 ). The vectors |0) and |0) are respectively the sl 2 highest and lowest weight states 
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and from (2.4) and (2.7) we readly obtain 



A(X,e)\0) = JJ/(A-^ + 7) |o> 

i=i 



(3.2) 



The framework. In order to explore the relations (3.1) and (3.2) we shall consider 
the quantity 



^(a ,# + 7) n + 1? - ^ 1°) 

i=i 



(3.3) 



computed in two different ways. We first compute (3.3) using (3.2) to find that it is 



proportional to Zg_ 7 (Ai, . . . , Al). Next we compute (3.3) using the relations (3.1) to 



move the operator A(X , 9 + 7) through the string of operators B(Xj, 8 + (j — 1)7) and 
then consider its action on the vector |0). By doing so we obtain a combination of terms 
Zg depending on the set of L + 1 variables {Aq, Ai, . . . , Xl} taken L at a time. 



Functional equation. The relations (3.2) follows from the slz highest weight repre- 



sentation theory and its consistency with the dynamical Yang-Baxter algebra implies the 
functional equation 



Mq Ze_ 7 (Ai, . . . , Xl) + Zg(Xo, . . ■ , Ai_i, Aj+i, . . . , Ax) 



i=0 



with coefficients given by 

L 



/(f) 



/(« + J>y) 



II f(Xo - ft) 



3=1 
f(0 + l) 



/(e + (L + i) 7 ) 



n/( A o-^+7)n 



j'=i 



/(Aj - A ) 



/ (g + 7 + Ao -A I ) /( 7 ) 
/(0 + (L + l) 7 ) /(Aj - A ) 



11/(^-^+7)11 



3=1 



/(A, - A, + 7) 
/(A, - A,) 



(3.4) 



(3.5) 

Some remarks are in order at this stage. Although the partition function considered 
here reduces to the one studied in [lj when the elliptic theta-function / degenerate into 
a trigonometric function, the functional equation (3.4) still differs significantly from the 



one obtained in |T|. For instance, ( |3.4[ ) is a functional equation also over the variable 9 
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and even in the limit 9 — > oo, where Za_~ and Zq coincide, 



26 



3.3) 



a functional equation different from the one presented in 
to the fact that here we have started our analysis with the quantity 
(0| C(X , 9 + 7) Uj=l B {\, + (j- 1)7) |0) as employed in the works |1 and 
different starting point allows us to obtain a simpler functional equation whose solution 
will be discussed in the next section. 



we still would be left with 
This divergence is due 
instead of 
This 
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4 The partition function 



This section is concerned with solving the functional relation (3.4). The method we 



shall employ is essentially the one described in [j] which exploits special zeroes of Zq 
to produce a separation of variables. Some structural properties of (3.4) will be of 



utility to help us identifying the elements required to solve this functional equation. For 
instance, the partition function Zq is a function of two sets of variables, i.e. {Ai, . . . , A^} 
and {/ti, . . . , /xz,}, in addition to the parameters 7, 9 and the elliptic nome p. In our 
framework, however, the set of variables {/xi, . . . , /xl} can also be regarded as parameters 
while 9 is promoted to a variable. This follows from the fact that (3.4) is an equation 
not only over variables \~ but also 9. 



With this in mind we can see that (3.4) is a homogeneous equation in the sense that 
aZg is a solution if so is Zq and a is independent of A,- and 9. This property implies 



that the equation (3.4) will be able to determine the partition function up to an overall 



multiplicative factor independent of Xj and 9 at most. Thus the complete determination 
of Zq will require that we are able to compute it for a particular value of Xj and 9 in 
order to determine this overall factor. Any point on the (Xj, 6 l )-plane would serve our 
need and we can choose it such that the evaluation of Zq is as simple as possible. As 
demonstrated in the Appendix |A| the evaluation of Zq in the limit (Xj,9) — > 00 can be 



performed in the same lines of 26 . Moreover, the equation (3.4) is linear which raises 



the issue of uniqueness of the solution since linear combinations of particular solutions 



also solves (3.4). Similarly to the case considered in we will see that the location of 
zeroes of Zq will select the appropriate solution uniquely. These properties are discussed 
in the Appendices A and B, and they are summarised as follows. 



Asymptotic behaviour. In the limit (Xj, 9) — > 00 the partition function (2.10) be- 
haves as 



Zg(Xi, . . . , Xl) ~ 

n / \ T OO OO 

f(l) L 



L (a) L(L-l) 



,,,-tj E E - E ■■■ E (-i)^"'" 

L L-l 



(i) (i) (-£-) (£) 

n\ ——00 n^_ 1 =— 00 — — oo n^J 1 =— oo 



nn 

a=l i=l 



(4.1; 



where e n = e ( 2n+1 \ p n = ^("+2) ; q n = and /i^ = : i 7^ a}. Here Sl denotes the 
group of permutations of L objects and a = o~(l) . . . cr(L) stands for a given permutation. 



7 



The set of inversion vertices for a given a is denoted by I a . 



Higher order theta-function. The partition function Zq is a theta function of order 
L and norm tj in each one of its variables Aj separately. That means it can be factorised 

as 



(0 



where C is \ independent and the zeroes £ w satisfy ^f 1 + 



(4.2) 

ti. Unveiling 



special zeroes ^f 1 will be an important step for solving (|3.4l). 



Now we proceed with the analysis of (3.4) in the lines of [T]. For that we look for 
special values of the variables Xj such that particular zeroes of Zq can be identified. 



Special zeroes. The coefficients M and given in (3.5) exhibit a factorised form 
and due to that identifying their zeroes is a simple task. For instance, when Ao = /xi and 
Ai = fix — 7 we find that Mo = Nq = N\ = 0. Next we can set Xj = Xj+i — 7 successively 
for j e [2,L — 1] and collect the result at each step. At the last step we find 

iV 2 Zffi/ix, n x - 7, X L - (L - 2)7, X L — (L — 3)7, . . . , X L ) 







(4.3) 



and since N 2 is finite we can conclude that the vanishing of (4.3) is due to Zq. This 



result can now be substituted back into the previous steps leading to (4.3). By doing 
so we find the more general vanishing condition, namely Zq(hi, /ji — 7, A 3 , . . . , A^,) = 0, 
for general values of the variables Xj with j e [3, L]. Now this process can be repeated 
starting with variables Ao = \i\ and Aj = \i\ — 7 for any j e which allows us to 

conclude that Zq{[1\, . . . , fix — 7, . . . ) = 0. 



Building up the solution. The zeroes of Zq above unveiled have a special appeal 



since we are interested in the solution of (3.4) consisting of a higher order theta-function 



(4.2). Taking that into account, those special zeroes imply that 



Z e (fii, A 2 , . . . , A L ) = Yl f(^j - A*i + 7) V e (X 2 , . ■ ■ , X L ) 
i=2 



(4.4) 



where Vq is also a theta-function but of order L — 1 in each one of its variables. Next 



we substitute the expression (4.4) back into the equation (3.4) and set A = /ii- The 



resulting equation can then be solved for Zq(Xi, . . . , Xl) yielding the expression 



Zq(Xi, . . . , Xl) — 2_j m * • • ' ^i-ij ^t+ij • • • 



i=l 



(4.5) 



S 



with coefficients 



rxii 



f(9 + 7 + nx - Xi) -A- f{Xi - iH + 7) -tt f( x 



f{o + i) 



3 . =2 /(Mi-Mi + 7) j=1 

7^ 



Ml) 



/(A, -A, + 7) 
/(A, - A,) 



(4.6) 



We then substitute the formula (4.5) back into the original equation (3.4) and set 



Hi. After eliminating an overall factor we are left with the equation 



L-l 



Pa Ve- y (Xi, . . . , Xl-i) + Qi Vg(Ao, . . . , A 



i— 1) Aj+i, 



X L -i) =0 , 



(4.7) 



8=0 



where the coefficients Pq and Qi correspond respectively to the coefficients M and iV, 
given in (3.5) under the mapping L — > L — 1, 9 — >■ + 7 and /ij — >■ //j+i. Thus the 



function Vg obeys essentially the same equation as the partition function Zg but for a 
square lattice of dimensions (L — 1) x (L — 1). Now since VJ? is also a theta-function, 
this procedure can be repeatedly carried out until we reach the equation for L — 1. The 



solution of (3.4) for L = 1 can be found in the Appendix D and gathering our results we 



obtain the following solution for general L, 

Zg{Xi, . . . , Xl) 

where 



P(r(i)...(T(L) 

&es L 



V(l)...<r(L) 



L f{9 + wy - A CT ( n) + fj^) 



nL? /(mi - ^ + 7) n _. 



n 



/(fl + n 7 ) 



n 



m>n 



f{X a (m) ~~ Acr(n) 



+ 7) 



/(A 



c(m) 



A CT („)) 



n /(A CT (n) - m? + 7) n /( a ^w - Mi) 

(4.9) 



j'<n 



The o vera ll factor Ql arises from the homogeneity of (3.4) as previously discussed, and 
from (4.1 ) we obtain Q L = f( r y) L Ylt=2 /(Mi — Mfc + 7)- ft is important to remark here 



that this partition function has also been considered in [16 18 where a similar but still 
different expression for F a (i)...a(L) has been found. 

Multiple integral formula. The partition function Zg can be represented by a mul- 



tiple contour integral as follows. The function Vg in the formula (4.5) is essentially our 
partition function for a lattice of size (L — 1) x (L — 1) and modified parameters. In 



fact, the decomposition of Zg in terms of Vg as described by (4.5) can be thought of as 



a separation of variables. Moreover, the prescription given by (4.5) can be mimicked by 
the Cauchy like integral 



Zg(Xi, . . . , Xl) 



H l ({X 3 }\{w 3 }) 

r£,-=i f(w* - a 3 -; 



yrdWj 



2ivr 



(4.10) 
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with integration contours enclosing solely the zeroes of / when W{ — > Xj. Also we assume 
that H L ({Xj}\{vjj}) has no poles inside the integration contour. 
For the case L = 1, we can immediately read from (D.6) that 



fMAxK) = /'(o)/(7) 



f(6 + 7 - wt + fix) 



(4.11) 



f(0 + l) 

where f'(0) denotes the derivative of f(X) with respect to A at the point A = 0. Guided 
by this result when looking to (4.5) we then find 

L L 

^({a^ik}) = [f(o)/(7)] L n - ^ + 7)/k - w t ) n 



fie+jj-Wj + iij) 



f(o + n) 



X 



n - w ^ n - Pi + 7 ) • 



(4.12) 



The expression (4.12) already takes into account the asymptotic behaviour (4.1) and 
though here we have considered a functional equation different from the one obtained 
in [l], the expression (4.12) indeed reduces to the formula of [T] in the degenerated limit. 
Moreover, it is worth remarking that the homogeneous limit Xj — > X and fij — > can be 
trivially obtained from the integral formula (4.10, 4.12). 



5 Concluding remarks 



In this work the partition function of the elliptic SOS model with domain wall boundaries 
was studied through a fusion of algebraic and functional techniques. The partition 
function of the model was shown to obey a functional equation provenient from the 
dynamical Yang-Baxter algebra which is valid for general values of the model parameters. 
The solution was then obtained as a multiple contour integral. 

The possibility of deriving functional equations for such partition functions from the 
Yang-Baxter algebra and its dynamical counterpart was firstly demonstrated in 26 ,1]. 
Although here we have also employed the dynamical Yang-Baxter algebra, the mechanism 
considered in Section [3] differs from the one used in 26 T] , and the resulting functional 
equation is significantly simpler than the ones previously obtained. Interestingly, solving 
this new type of functional equation follows the same lines of fl] but each one of the 
steps required are dramatically simplified. 

The elliptic SOS model considered here is also referred to as 8VSOS model in the 
literature and for the special value of the anisotropy parameter 'y = it reduces to 
the so called Three- colouring model 27 . For the case with domain wall boundaries, the 
partition function of the Three-colouring model was shown to obey a certain functional 

In 



equation in 
the work 



21 



28 , 29 but a possible connection with our results has eluded us so far. 



this same partition function was studied under the light of the symmetric 
theory where a set of two- variables polynomials have been introduced. These 

to satisfy a certain partial differential equation and 



polynomials 

polynomials were conjectured in 



30 



recurrence relation, to which (3.4), (4.10) and (4.12) might shed some light into their 
proofs. 
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A Asymptotic behaviour 



In the limit 9 — > oo the TZ- matrix (2.4) resembles the one associated with the six- vertex 



model, except that the Boltzmann weights (2.5) still consist of elliptic theta-functions. 



Moreover, from the definition (2.10) we can readly see that the whole dependence of Zq 
with a particular variable Xj will be described by the operator B(\j,6 which is 

very similar to the six- vertex model analogous in the proposed limit. Thus considering 



the definition (2.3), in the limit (Xj,0) — > oo we obtain 



B{\ v e + 31 ) 



2 L- 



L-l 



ni=— oo ni,_i=- oo 



(L-i) 
2 



i=l 



c ni • • • c nj_i 1 j z-rij ■ ■ ■ c n L _! ■ 



(a.i; 



In the expression (A.I) we have introduced the conventions e n = e ( 2n+1 ) ; p n = p 



(n+i) 2 



q n = el and n 



[ni, 



^l-i), while the operator P n reads 



P n = K 



K' 1 



(A.2) 



In their turn the operators K n and X are defined as 



ql 
q n 



X~ 





1 



(A.3) 



and, together with X + 
SU-2 algebra 



1 




, they satisfy the following analogous of the g-deformed 



±70—1 



K n X K, 



q n + m X 



± 



2 

K n K m - K~ x Kz} 



n m 



(q n+m — q 



-1 - 

n-\-m 
2 



(A.4) 



The relations (A.4) allow us to demonstrate the properties 



pn("> pnW 
i 3 

i i 







(» < 3) 



(A.5) 
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which will be of utility for the analysis of the asymptotic behaviour of Zg. 

In order to analyse the behaviour of Zq in the proposed limit we now substitute the 



expansion (A.l) into the definition (2.10), and use the relations (A. 5) to reorganise the 



result properly. By doing so we obtain 

Zg(\i, . . . , X L ) ~ 
f(l) L 



2HL-1) 



it E - E ■■■ E ■■■ E 



-l~)Ea=iE l= i n 



L-l (a) L(L-l) 



-oo n 



(l) 



(i) 

oo n, — — oo n 



(L) 



*L-\- 

a=l i=l ctScSl (a,b)elcr a=1 

(A.6) 

where /i^ a ^ = {//j : i ^ a}. As usual Sl denotes the group of permutations of L objects 
and a(a) stands for the permutation of the a-th object. In order to clarify the meaning 
of l a let us consider the usual two row representation of a. We draw a line starting at 
the object a in the top row and ending in the bottom row at the position a (a) such that 
only two lines intersect at any one point. The lines are labelled by their numbers in 
the top row and the points of intersection are called inversion vertices. In this way an 
inversion vertex can be labelled by a pair (a, b) with a < b such that a and b label the 
two intersecting lines originating the inversion vertex. Then denoting [L\ = {1, . . . , L}, 
we call I a = {(a, b) G [L] x [L] : a < b and a(a) > <j(b)} the set of inversion vertices 
labels of a given permutation a. 



The next step to obtain an explicit expression for (A.6) is to compute the quantity 



Tla=i -P^ (a) |0) which can be readly performed since the operators P™ w consist of 



(a) 



simple tensor product (A. 2). Thus considering (2.11) we obtain 



L L-l 



<oinc w io>=nn4> 



(A.7) 



0=1 



a=l i=l 



which can be substituted in (A.6) yielding 



Z$(Xx, . . . , Xl) 
f(l) L v 



E •■■ E 



E (-i) 2 -' 2 



£ ~1 „(") L(L-l) 



(1) (1) 
n l =— OO n L — i 



L L-l 



(a) 



a=l i=l 



e n 1 

a£S L (a,b)el„ 



(A. 



in the limit (A^, 6) — > oo. 
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B Theta-function properties 



In this appendix we recall some useful properties of elliptic theta-functions that we have 
considered through this paper. We remark here that many of these properties have also 
been discussed in [16]. The function / defined in Section [2] consists basically of the 
Jacobi theta-function Q\ 25 and in this paper we have ommited the dependence of / 



with the elliptic nome p for brevity. In what follows we summarise some properties of 
elliptic theta-functions adjusted to our conventions. 

Addition rule. The function / satisfy the addition rule 
/(Ai + A 2 )/(A! - A 2 )/(A 3 + A 4 )/(A 3 - A 4 ) = 

/(Ai + A 4 )/(A! - A 4 )/(A 3 + A 2 )/(A 3 - A 2 ) + f{X, + A 3 )/(A 1 - A 3 )/(A 2 + A 4 )/(A 2 - A 4 ) . 

(B.l) 



Analiticity and periodicity. The function / is an entire function, that is to say all 
of its singularities are removable, and it has only simple zeroes. It is also an odd function 
and quasi doubly-periodic, i.e. 



/(A - hr) = -/(A) 



/(A - ivrr) = -e 



2A — YKT 



/(A) 



(B.2) 



Trigonometric limit. In the limit p — > we observe the degeneracy lim p ^.o/(A) 



sinh(A). This fact motivates the use of the definition (2.3) which allows for an easy 
comparison with previous results in the literature. 



Higher order theta-functions. For a fixed value of the elliptic nome r, we call J 7 a 
theta-function of order L and norm t if 



(B.3) 



for constants Xj such that Yuj=i Xj — t- Moreover, due to (B.2) one can readly show the 
quasi-periodicity 



F{\ - ivr) 
F(\ - iTrr) 



■l)^(A) 

■l) L e 2(LA - t) - i7rr - L J"(A) 



(B.4) 



In fact, the factorised form (B.3) and the quasi-periodicity ( B.4[ ) for entire functions can 
be shown to be equivalent properties 31 . This feature allows us to state a more general 
result. Let T be defined as 



t j=i 



(0^ 

3 > 



(B.5) 



such that ^2^ =1 Xj = t f° r an y The function T is entire and obeys the quasi-periodicity 



(B.5), thus it can be factorised similarly to (B.3). 
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Zq as a higher order theta-function. The partition function Zg defined in (2.10) 
is written as a product of operators B(X,9). As a matter of fact, the whole dependence 
of Zg with a particular variable Xj is contained in a single operator B(Xj,9 + jj) since 
the vectors |0) and |0) are constants. From the definitions (2.4), (2.5), (2.7) and (2.9) 
we can clearly see that the entries of B(Xj, 9 + jj) will be of the form (B.5). Thus the 
partition function Zg is a theta-function in the variable Xj of order L and norm tj. 



C Zq as a symmetric function 



The commutativity of operators B(X, 9) as described by (3.1 ), together with the definition 
(2.10), implies that Zq is a symmetric function. This commutativity has been extensively 
employed in the derivation of (3.4) and here we intend to show that this symmetry 
becomes an inherent property of the solutions of (3.4). 

Through the inspection of the coefficients (3.5), we notice that Ni <H- Nj under the 
mapping Xi -H- Xj while Mq — > Mo and Nk — >■ Nk for k ^ Thus performing this 
mapping on (3.4) and subtracting it from the original equation we obtain the following 
relation, 

Mo[Z0_j(Xi, . . . , Xi, . . . , Xj, . . . , Xl) — Zq_ 1 (Xx, . . . , Xj, . . . , Aj, . . . , Xl)\ 
+ Nq[Zq{Xi, . . . , Xi, . . . , Xj, . . . , Xl) — Zg{Xi, . . . , Xj , . . . , Xi, . . . , Xl)] 



+ NkZg(Xo, . . . , Xi, . . . , Afe_i, Xk+i, ■ ■ ■ , Xj, . . . 



— Nf.Zg(Xo, . . . , Aj, . . . , Afc_i, Afc + i, . . . , Xi, . . . ) — . 



(c.i; 



)c=i 



Next we solve (C.I) for the Z-th term of the summation over the index k which yields 
the expression 



Mo 

Nn 



[Zg- 1 (Xi, . . . , Xi, . . . , Xj, . . . , Xl) — Zg_^(Xi, . . . , Xj, . . . , Xi, . . . , Xl)} 



+ -irr[Ze(Xi, . . . , Aj, . . . , Aj, . . . , A^) — Z s (Xi, . . . , Xj, . . . , A,, . . . , X L )} 
+ ^ — ^6»(Ao, • • • , Xi, . . . , Xk-i, Afc+i, . . . , Xj, . . . ) 



k=l 



-irrZg(Xo, . . . , Xj, . . . , Xk-i, Xk+i, ■ ■ ■ , Xi, . . . ) 



fc=i 



Zb{Xq, . . . , Xj, . . . , Xi-i, Aj+i, . . . , Xi, . . . ) — Ze(Ao, . . . , Xi, . . . , A;_i, A; + i, . . . , Xj, . . . ) . 

(C.2) 
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The RHS of (C.2) does not depend on A; so this variable can be chosen such that the 

an c 

Aj, 



LHS of (C.2) vanishes. Thus we can conclude that 
Ze(Xo, . . . ,Xj, . . . , Xi-x, A^+i, 



> Aj, 



Az-i, A; + i, . . . , Aj, . . . ) , 
(C.3) 

and since this is valid for any i, j and I, the symmetry property 



Z e (. . . , Aj, . . . , \j, . . . ) — Zq(. . . , \j, . . . , Aj, . . . ) 



(C.4) 



immediately follows. 



D Solution for L = 1 



The functional equation (3.4) for L = 1 explicitly reads 



with coefficients 



M 



MoZ^Ax) + iVo^(Ai) + iVi^(Ao) 
f(0) 



(D.l) 



/(A -f*i) 



f{0 + l) 

f(e + 7) 

/(# + 7 + A - Ai 



/(A -/ii +7) 



/(7) 



/(0 + 2 7 ) /(A!-A ; 



/(A1-A0 + 7) 
/(Ax - A ) 

■/(Ai + 7) 



(D.2) 



By setting Ao = Ai — 9 — 7, the coefficient N± vanishes and (D.l) simplifies to 

f(0 + j) 



/(0 + 7 + ^-A, 



Zo-i(X\ 



(D.3) 



The relation (D.3) is an equation only over the variable # which is readly solved by 

/(0 + 7-Ai + pO 



^(Ai 



-F(A0 



(D.4) 



/(* + 7) 

where F is ^ independent. After eliminating the dependence with 0, we can substitute 



(D.4) back into (D.l). The resulting equation can then be simplified and we obtain the 



relation 

f(l)f(0 + J + m- AoXffj + 7 + Ap - Xi)f(\i - /ii + 7) 
/(0 + 7)/(0 + 2 7 )/(A 1 -A o ) 



( J P(Ai)-F(A o ))=0. (D.5) 



behaviour (A. 8). Thus we find for L = 1, 



From (]D.5|) we can conclude that F is a constant and it can be fixed by the asymptotic 

find for L = 1, 

MX) = m f{6+ Z- X + ^ . (D.6) 
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